Spin-| Collective Excitations in BEC of Interacting Spin-1 Atoms 
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We construct spin-i collective excitations in BEC of interacting spin-1 atoms. These excitations 
exist in states with a maximal global degeneracy. The stability and energy of these objects are 
determined by interactions with spin fluctuations and are studied based on a duality relation between 
hyper-monopoles and magnetic monopoles in physical space. 



Theories for Bose-Einstein condensates(BEC) of inter- 
acting spin-1 atoms were recently investigated by a few 
groups motivated by experiments on spinor BEC 

. Of particular interest is spinor BEC with two-body 
scatterings dominated by the total spin equal to 2 chan- 
nel, i.e., with an antiferromagnetic interaction. Com- 
pared with BEC of spinless atoms, this problem poses 
greater challenge to theorists and attracts considerable 
interest because of the following aspects. First, the 
ground state of noninteracting spin-1 bosons, once con- 
densed, has a spin degeneracy proportional to the square 
of the number of bosons involved. A spin dependent two- 
body interaction should be taken into account for the un- 
derstanding of the low energy properties Second, 
a spin long range order could be present in the thermal 
dynamical limit. Scattering between Goldstone modes, 
however, disrupts the order in the system and is crucial 
to the stability of the rotation symmetry broken BEC 
state. Third feature here is that the order parameter 
space IS [U(l) X S^]/Z2 @. This entanglement of a unit 
circle and a unit sphere was shown to lead to Z2 gauge 
fields and the existence of other exotic many-body states 
such as a condensate of singlet pairs 

Though a qualitative understanding of the nature of 
many-body ground states in spinor BEC with antiferro- 
magnetic interactions has been achieved after these the- 
oretical efforts, the nature of collective spin excitations 
hasn't been explored thoroughly and remains to be un- 
derstood. The purpose of this paper is to illustrate the 
possibility of having spin-i excitations, with their quan- 
tum number being fractions of that of an atom in BEC. 
Main results of this paper are that 1) a spin-i collective 
excitation carries a hyper-monopole defined in an inter- 
nal space, which can be seen as a magnetic monopole 
with respect to a hedgehog defect in spinor BEC; 2)the 
energy of a spin-i excitation is determined by interac- 
tions between hypermonopoles and spin fluctuations and 
can be finite in the thermal dynamical limit. 

To demonstrate the existence of spin-i excitations, we 
are going to introduce a lattice model. At each site, a 
spinor BEC of A^o spin-1 atoms with antiferromagnetic 
interaction is formed and is described by a Hamiltonian 
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ao)/3M; 00,2 are the scattering lengths of two bosons in 
total spin F = 0, 2 channel respectively; M is the mass 
of an atom. A'o(>> 1) is the number of bosons at chem- 
ical potential /i = pgjy. p is the number density of atoms; 
A^ and L are the number and total spin operators of the 
spinor BEC. Eq.l is valid as far as the spinor BEC can 
be treated as a rigid rotor. 

As shown in previous papers L can be consid- 

ered as an angular momentum operator defined on a unit 
sphere hyperspace of n, a nematic director. L and 
n satisfy a commutation relation: [Lcii^] = iheaP'j'n.'y 
(ea/37 is an antisymmetric tensor). Similarly, A^ is an an- 
gular momentum operator defined in a hyperspace, a unit 
circle exp(ix) and [A'^, x] = i^- In a representation of |n > 
and Ix >, eigenstates of collective operators n and Xj the 
collective excitations are characterized by a wavefunction 
^'(n, x). The symmetry of the N-boson wavefunction re- 
quires n,x + tt) = \E'(n, x). If an arbitrary state is 
labeled with two quantum numbers of operators A^ — Aq 
and L as (n, I), a physical excitation has to be of a form 
(n, 2/ — n), with n, / as integers. 

Most generally, the Hilbert space defined at each site 
by the Hamlitonian in Eq.l consists topologically dis- 
tinct sectors. To classify them, we then introduce the 
following global gauge transformation |n >^ g;i(n)|n >, 
Ix >^ 9v{x)\x >■ The first one defines an fiber bun- 
dle over a unit sphere 5'^ (the base space), i.e., E : S'^xS'^ 
specified by a mapping gh{n) : ^ S"^. The second one 
is for a usual U{1) bundle over a unit circle specified by a 
mapping gv{x) ■ ^ S^. In both cases, gauge fields can 
be topological nontrivial. Namely they are characterized 
by hyper-monopoles and hyper-fiux lines with charges 
given by winding numbers qh,v 
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Here gs = 47r?i^(a2 - ao)/3M (> 0) and gp = 'iiih^{2a2 - 



here 6*1^2 — and n = (sin cos 0, sin sin ip, cos 9). Tr 
is carried over a spinor representation of gh{n) (specified 
after Eq.3.). Each topologically distinct sector defined in 
Eq.2 is characterized by qh^qv, (lh,v = 0,±1,±2,.... 

L and A^ can be expressed as differential operators de- 
fined in hyperspace 5*^ and S^. Particularly, in sector 
qh , qv , the gauge invariant forms of L and A^ are 
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L = nx[*a^-A'^(n)],iV = ^--A'^(x) (3) 

where we introduce gauge potentials A''(n) = 
iTrg~^{n)dngh{n), A" = ig'^d^g^ in hyperspaces 5*^ 
and 5*^ respectively. Within each sector, the gauge po- 
tentials differ by a pure gauge transformation which pre- 
serves the invariant defined in Eq.2. For sector qh — 1, or 
more conveniently a hypermonopole of unit charge, the 
singular gauge transformation g/!.(n) — u{n) is defined 
as a spinor in a 2d representation: (cos6'/2exp(— j(/)/2), 

sin (?/2 exp(i0/2)) , corresponding to have two half Dirac 
strings inserted at the south and north pole of the unit 
sphere. The vector potential on the unit sphere A'' — 
|(n X e^)(n • ez)/\n x e^p. For a sector g„ = 1, or a 
hypcrflux of unit charge, the gauge potentials arc A" = 1 
and A'* = 0. 

Following Eqs.1-3, it is clear that in the absence of a 
hypermonopole, one can choose A'* — and redrive the 
excitation spectrum obtained in previous work. However, 
in sector = 1 where a hypermonopole is inserted, the 
angular momentum is quantized at half integers as in the 
Dirac string problem Spin-^ excitations emerge as a 
result of the insertion and their wavcfunctions are 

^'_i(n) = sin^exp[-i^],*i(n) = cos^exp[i^]. (4) 

These exotic excitations are surely forbidden when 
each site is decoupled from each other, following the con- 
sideration of the symmetry of N-boson wavcfunctions at 
each site. An individual hypermonopole does not exist 
in this limit; say it differently, it is always paired with an 
anti-hypermonopole to form a bound state within each 
site. However, when coupling between the BEC at differ- 
ent sites is finite, the symmetry of bosonic wavefunction 
should be imposed only globally instead of locally at each 
site. Meanwhile, the total spin is a conserved quantum 
number and sectors such as qh — I, which are eigenstates 
of the total spin, should be taken into account for the 
study of spin excitations. Thus, it becomes possible to 
have a hypermonopole localized at one site (or a region of 
a few sites) which is spatially well separated from the an- 
tihypermonopole at the other site. In a region where the 
hypermonopole is pinned, a collective excitation earring 
precisely one half spin appears. The interaction between 
spin-i excitations is mediated by spin fluctuations. We 
should prove in a 3d lattice, the interaction between these 
spin-i excitations can be weak and these collective exci- 
tations are well deflned objects. 

The estimate of the stability of a spin-i excitation and 
the energy is based on a duality between hypermonopoles 
and "magnetic" monopoles in real space. It is convenient 
to introduce a connection field A(r) = e^^e^Ve^ in real 
space, where and n form a local orthogonal basis; 
Gq • = Saf3 and • n = 0. The field strength is a 
gauge invariant Pontryagin density Hfc = e^^'^djAi = e^^^ 



gaPTji^ d^n^djii^. A magnetic monopole carrying charge 
Qm located at R™ is defined as a topologically distinct 
configuration of n(r) where 

Q,„ = ^ dr • H, or V-H = Q„,5{y - R"); (5) 

Qm ~ ±1,±2,.... All configurations with a monopole 
Qm located at R™ differ only by a gauge transformation 
in the connection field A(r) A(r) -I- f;~^(r)Vgm(r), 
gm(j) is an analytical unitary function. 

$(r"; R"*) as the many-body wavefunction in the pres- 
ence of a monopole at site R™ and a spin-i excitation 
or a hypermonopole at site r" can be written as 

<i>^r *o({n(r^)})vI/i(n(r"))|{n(r^)}) 

No 

|{n(r'')})^.9,(n(r")) J] l[ui{n{r-'))vi{n{rr)) 

7=/3,a 1 = 1 

(6) 

The summation is carried over all configurations in sec- 
tor Qm = 1 with monopole located at R'". ^'o({n(r''')}) 
is an arbitrary single valued function with 7 labeling all 
sites including a; (3 labels all lattice sites expcept site a. 
We have introduced P as a permutation operator to pre- 
serve the symmetry of the wavefunction. In a standard 
three-component representation of F = 1 state, uv{n) 

= ( — sin6'exp(— z(/))/-\/2, cos 6*, sin6'exp(i(/))/\/2)"'". 

Let us first consider a state where \E'o({n(r''')}) = 
(5(n(r'^) - n™(r'' - R")). Here, n™(r) = r/|r|. Follow- 
ing Eq.5, this hedgehog carries a monopole field H™(r) = 
r/|rp and satisfies Eq.5 at Qm — 1. By examing Eq.6, 
we prove that the Berry's phase Tm{h){C) of moving a 
hedgehog at r" (a hypermonopole at R™) around a hy- 
permonopole (a hedgehog) along loop C is 

r™(C)--r;,(C)[mod27r] = ^^. (7) 

il{C) is the solid angle spanned by loop C with respect 
to the hypermonopole (the hedgehog). Eq.7 shows that 
hedgehogs and hypermonopoles see each other as "mag- 
netic" monopoles with one half monopole strength but of 
opposite signs. 

One should notice that the Berry's phase of mov- 
ing a hypermonopole (a hedgehog) along path C van- 
ishes identically in the absence of a hedgehog (a hyper- 
monopole). It also follows that the wavefunction of two 
hypermonopoles at r"^ and r"^ is symmetric with respect 
to interchange of positions, 

$(r"Sr"=) = $(r"%r"i). (8) 

The exotic excitation therefore carries spin-^ but obeys 
bosonic statistics 0. 
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In a 2d lattice, a Skyrmion texture can be defined 
by choosing vE'o({n(r''')}) = S(^n{r'') — n''(r''' — R*)) in 
Eq.6. For a skymion located at R^ with n*(r) = 
(sin 6'(|r|) cos 0(r), sin 6'(|r|) sin 0(r), cos 6'(|r|)) and a hy- 
permonopole located at r", the Berry's phase Th of mov- 
ing a hypermonopole around a loop C of radius p centered 
at Rs is again connected with Ts of moving a Skyrmion 
around a hypermonopole along the same path, 

r,(C) = -r,,(C)[mod27r] ^ tt{1 - cos0{p)) . (9) 

We have assumed (j>{r) = argVy/rx, and 0{\r\) is an ar- 
bitary function whose asymptotic behaviors are given as 
6l(|r| = 0) = and 6l(|r| >> Co) = tt. The Berry's phase 
Ts in Eq.9 varies from at p = to 27r at p = +oo | |lO[ |. 

Eqs. 7,9 are valid for any state ^o({n(r''')}) in Eq.6 
therefore in a disordered limit, if Tm,s are defined as the 
Berry's phases due to the presence of a hypermonopole. 
In terms of A(r) fields, Th in Eq.7 indicates that a 
hypermonopole with velocity experiences a Lorentz 
force iv/i X H™(r" — R™(t)) in the connection field of 
a hedgehog located at R™. Tm in Eq.7 on the other 
hand shows that a hedgehog with a velocity Vm experi- 
ences a Lorentz force iv„i x H'"(R™(i) — r") at R™, 
with H™ being a monopole field of a hypermonopole at 
r". So for the static hypermonopole, an electric field 
E(r") = -i(9fA(r" - R™(t)) is generated due to the 
motion of a hedgehog; A is precisely the connection po- 
tential of a magnetic monopole field H™ = V x A. We 
conclude a hypermonopole carries an electric charge with 
respect to potentials (Ai,^o) and is minimally coupled 
to the connection fields: ji • A^ -I- pAo. 

At this point, we should emphasis that with respect 
to hypermonopoles, connection fields of Q„i = ±1 are 
physically distinguishable, disregarding the Z2 symmetry 
found in previous works That is Hfc(n) = — Hfc(n), 

upon an inversion n — n and H((3„i = +1) = 
— H(Q„j = — 1) It basically refiects a conversion 

of a hypermonopole qh = I into an antihypermonopole 

= —1 upon the inversion of n, following Eq.6. 

However, a homotopical distinction between the con- 
nection fields of Qm = ±1 exists only when the influence 
of a linear defect is absent. Otherwise, a hedgehog with 
Qm = +1 transforms into a Qm = — 1 one when mov- 
ing around a 7r-disclination and half-vortex composite, 
or a Z2 string Moreover, a hypermonopole develops 
a Berry's phase Fs in a matrix form or Fs = <Tx7r/2, 
with ax an x-component of Pauli matrix acting on a 2d 
hypermonopole-anti-hypermonopole space. By moving a 
Qh = ^,Qm = 1 compound shown in Eq.6 around a Z2 
string, one obtains a, = —1, Q„i = 1 compound. So in 
the presence of Z2 strings, hypermonopoles of = ±1, 
or the connection field Hk of Qm = ±1 hedgehogs be- 
comes homotopically indistinguishable. We will neglect 
this infiuence of Z2 strings in the discussion of the sta- 
bility of spin-i excitations, for Z2 strings are gapped in 



the limit that interests us. As far as the Z2 symmetry 
is broken, the connection field of Qm = ±1 should be 
considered to be physically and homotopically distinct. 

The low energy dynamics of the Pontryagin field 
or magnetic monopole {R™} can be obtained by study- 
ing the effective theory of n. As shown in |^,^, the low 
energy spin dynamics in the spinor BEG is described by 
a nonlinear sigma model Cs = 1/2/(9^^11)^. Derivative 
d^, p = 0,1,2,3 is defined in d + 1 Euclidean space 
= (f , f ); f = r/a, and f = itcs/a, with a as the lattice 
constant and Cg is the spin-wave velocity. In the current 
situation, 1/2/ = Et / Eso and Et is an exchange in- 
tegral between two adjacent sites, Eso — P9s/2No is the 
zero point rotation energy; and Cg = a\/EsoET- In spin 
disordered limit, the effective theory for the gauge field 
F/ii^(r) = e°'^''€f_,i,nadf_in/3di,nj can be obtained upon in- 
tegration of short range spin fiuctuations and was shown 
to have a form Cg = \/2g¥ , g ^ As being the 
spin gap in the excitation spectrum measured in units of 
Eso 10 • We will only be interested in a phase coherent 
but spin disordered BEG, where the connection fields are 
induced by spin fiuctuations as discussed above. 

A static spin-i excitation, which is represented by 
^'i(n(r"))|0 > and carries a hypermonopole, inter- 
acts with the gauge field via Hj — ^ J drS{r — 
r")Ao(r). The energy of a spin-i excitation therefore 

is = -Lo^^ ln(rexp (-i/g^" Hi{T)df)), with Lq the 
perimeter along temporal direction and being infinity. In 
terms of induced connection fields, 

E = -^\n{exp{-l [ dxo [ dxzFoa)).^ (10) 
^0 ^ Jo Jo ^'f-'-r 

and the average is taken over the partition function 
Z{{F^u}). Only the y-component of the electric field 
defined by F02 contributes to the energy. The calcula- 
tion can be carried out parallel to that in |l^ and we 
summarize results here. 

In {2 + l)d, merons are condensed; the connection field 
becomes massive and has short range correlations. The 
expectation value in Eq.lO, which is equal to the Wilson 
loop integral of a gauge field, decays exponentially with 
an exponent proportional to PLqL2', P = exjpl—Tr/g) is 
the quantum tunneling amplitude of merons. The energy 
of a hypermonopole is linear in term of the size of system 
in 2d, or spin-^ excitations interact via a linear potential. 

In (3 + l)d case, however, magnetic monopoles are 
gapped in weakly disorder limit (5 <C 1) because P 
vanishes as Lq goes to infinity. The connection field 
is massless and has a quasi-long range correlation. Hy- 
permonopoles or spin- i excitations interact with each 
other only through massless photons and spin-^ objects 
in BEG (spin disordered but phase remains coherent) are 
well defined free excitations. The energy of these exci- 
tations is comparable to the spin gap. In strongly dis- 
ordered limit ( g ^ 1), the action is finite and magnetic 
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monopoles again become condensed; situation becomes 
similar to that in (2 + l)d. What happens here is that as 
an electric charge, or spin-i excitation is inserted into a 
magnetic monopole condensate, the electric flux emitted 
forms a flux tube. In strongly disordered limit, the inter- 
action between two spin-i excitations is linearly propor- 
tional to the length of the electric flux tube connecting 
them. Only when monopoles uncondense as in a weakly 
disordered limit in a 3d lattice, an individual spin-i col- 
lective excitation is liberated from its interaction with 
the connection field. 

The condensation of monopoles also implies a reduc- 
tion of the ground state degeneracy. Monopoles inter- 
act with each other via a Coulomb force and form a 
Coulomb gas once condensed. In a Sine-Gordon rep- 
resentation for the monopole gas, or a Coulomb gas, 
where the monopole charge Qm couples with the Sine- 
Gordon field Xsg with a term exp{ixsgQm)i the parti- 
tion function Z = J Dxsg exp(— S'sg) and the action is 
SsG ~ — / dfP cos Xsg in a zero mode limit . Clearly, 
at the mean field approximation, the ground state degen- 
eracy is infinite when monopoles are gapped (P = 0) and 
singlefold (up to a 27r phase shift) when they condense 
(P 7^ 0) In addition, all BEG states discussed in 

this paper have twofold Z2 degeneracy in a cylindrical 
geometry, with .Z2-vortices gapped HJl^. Apparently, 
spin-i excitations exist only in states with a maximal 
global degeneracy. 

To conclude, we want to make the following remarks 
concerning the topological spin-^ excitations. First, spin- 
i excitation exists when spin is gapped and incompress- 
ible. For BEC of a finite size, the spin-^ excitation should 
be understood as a point like object localized at a hy- 
permonopole plus a uniform spin density along the edge 
of the BEC so that total spin of the system is an in- 
teger. This is similar to the charge fractionalization in 
fractional quantum Hall system [Q. Second, through a 
Hopf projection from CV^ complex fields of z = (zi,Z2) 
(z+z = 1) living on to a unit vector n on S*^, the 
NLcrM derived for spinor BEC in d can be formulated 
in a complex field representation. Spin-^ z-quanta inter- 
act with each other through compact C/(l) gauge fields 
which have confining-deconfining phases in (3 + [ p^ . 
The explicit construction of spin-i excitation shown here 
suggests z-quanta as hypermonopoles. Third, spin-i ex- 
citations are also believed to present in an antifcrromag- 
netic spin lattice with integer underlying spins. For an 
5=1 antiferromagnetic spin chain, with Hamiltonian 
H — ^ .[Si-Si+i— /^(Si-Si+i)^], a valence-bond sofid with 
onefold degeneracy is expected to be degenerate with a 
dimerized state of twofold degeneracy aX [3 — \ point 
I psf . The Bethe Ansatz solution (3 — 1 point reveals 
that one-particle excitations in the spin chain are spin-^ 



for spin-i antiferromagnetic systems [Q. In 2D frus- 
trated antiferromagnets, due to a Higgs mechanism, free 
spin-i spinous are argued to exist disregarding the un- 
derlying spins at lattice sites The hypermonopoles 
discussed here appear to be intimately connected with 
those fractionalized objects discovered in low dimension 
integer-spin antiferromagnetic lattices. 

Finally, a similar discussion for (1, 0) or (0, 1) collective 
excitations can be carried out. The existence of particle 
(1,0), which carries = 1 unit hyperfiux depends on 
whether the Z2 gauge field is confining or dcconfining, 
as suggested in I would like to thank F.A. Bais, E. 
Demler, F. D. Haldane, G. 't Hooft, R. Moessner, M. 
Parikh, J. Strict for discussions, particularly S. Kivelson 
for sending me his unpublished paper, and T. Senthil for 
useful comments on the manuscript. 
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